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ABSTRACT 


A comprehensive kinetic model is presented which accounts for 
the major and side reactions present in polyaraic acid and 
polyimide reactors. The mathematical framework is quite general, 
and permits incorporation of heat and mass transfer aspects 
typically encountered in large scale reactors. A few rate 
constants have been curve- fitted using some (limited) experimental 
data available in the open literature. Using these (and the other 
rate constants from the literature) , some interesting phenomena 
have been predicted at higher temperatures. Detailed experimental 
data are required to evaluated all the rate constants with 
precision . 
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CHAPTER 1 


INTRODUCTION 

Polyimides are an important class of speciality polymers, 
exhibiting several desirable thermal and mechanical properties^ 
These include excellent electrical insulating capability and 
unusually high thermo -oxidative stability. Polyimides can be 
prepared by a number of reaction routes . The most common among 
these is by the reaction between tetrabasic acid dianhydrides 

(e.g., pyromellitic dianhydride) with diamines (e.g., 

oxydianiline) in a two step process. In the first step, the two 

monomers react with each other in a suitable solvent [e.g., 
dimethylf ormamide (DMF) ] at about 25°C to give polyamic acid (PAA) 
precursors. Films or moldings of PAA are then cured at relatively 
higher temperatures of about 120-130°C. During this second stage, 
the PAA is converted to the polyimide by the dehydrocyclization 

reaction, while simultaneously releasing water as a byproduct. 

7 

since polyimides are relatively nonporous materials , this could 
lead to several technological problems like diffusional 

limitations (at a 'micro' level) of the reaction rate, crack and 
void formation, etc. The rate of volatile formation by chemical 
reaction must be synchronized with its rate of removal by mass 
transfer (at the 'macro' level), so as to minimize these problems. 
In addition, the water formed hydrolyzes the PAA molecules present 
and lowers the average molecular weight of the product . The 
properties of the polyimide product are, thus, affected by the 
average molecular weight of the PAA precursor, as well as the rate 
of the hydrolysis (and other) reactions. A proper design of these 

.*11 



reactors requires an understanding and modeling of the various 
physico-chemical phenomena associated with both these steps of 
polyimide manufacture. This study is an attempt along this 
direction. 

An improved methodology of studying this system is discussed 

here, which is based almost totally on the kinetic approach. A 

kinetic scheme is presented (see Table 1) , which incorporates the 

major reactions, as well as the important side reactions taking 

place in both the steps of polyimide manufacture. Mass balance and 

moment equations are written and appropriate moment-closure 
8 9 

conditions ' are developed. These equations are solved for an 

isothermal, well mixed batch reactor (first stage reactor in which 

the precursor is formed primarily) . The Box^*^ complex technique is 

then used to 'tune' those rate constants to which model results 

are most sensitive. This is done to minimize the sum of square of 

errors between the model predictions and experimental results of 
11 

Kolegov et al . for PAA synthesis . 

The kinetic model developed herein provides a general 
mathematical framework for both the first and second (polyimide 
formation) stages of reaction, even though the model is tuned 
using data for the first stage only. Qualitative trends for the 
second stage (in absence of mass transfer limitations) have been 
predicted to illustrate the general nature of the model and its 
applicability at higher temperatures. Additional tuning of the 
other parameters against data may need to be performed once 
experimental results for the second stage become available in the 
open literature. 

It should be added that this kinetic model can easily be 
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extended to predict some very interesting experimental 
observations described (qualitatively) by Johnson"^ such as the 
dependence of the average chain length of the product on the way 
in which the two solid monomers and the solvent are mixed in the 
reactor [e.g., higher molecular weight products are formed when 
the two solid monomers are blended in stoichiometric proportions 
and then added to a tank of pure solvent, than when either monomer 
(solid) is added slowly to a solution of the other monomer in 
solvent, etc.] . Such operations (involving mass transfer 
limitations) cannot be as easily modeled using the currently 
available theories , which involve varying amounts of 
probabilistic concepts. 
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CHAPTER 2 


FORMULATION 

The thermodynamics and kinetics of polyimide preparation have 

been extensively studied."^' Table 1 gives the kinetic 

scheme incorporating most of the important reactions occurring 

during polymerization. This scheme is more general than that used 

by Kolegov et al. (who omit the reactions asssociated with rate 

constants k^ , k^ , k^ and k^ in their kinetic scheme) , but is the 

same as the recent scheme used by Tsiang^^ (who use a different 
22 23 

approach ' at modeling the system) . The various functional 

groups, -A, -B and -A', at the ends of the molecules and internal 

* # 

linkages, -AB-, -A B- and -A B-, are also defined in this table. 
Since PAA can be present in two kinetically non-equivalent 
forms^^'^^' a transformation from one to the other (Rxn. 5) 

is incorporated in this scheme. Only the -AB- internal linkage is 
able to cyclize in the curing stage through Rxn. 2. Rxns . 2-5 

would be significant primarily in the curing stage, while Rxn. 1 
would be most important in the first stage. The kinetic scheme is, 
thus, observed to include all the reactions (relevant to both 
stages) . 

The rate constants used in Table 1 are given by the Arrhenius 
equation 

k. (or kM= A.exp(-E. /RT) (l) 

1 1 1 1 

The (guess) values of the frequency factors, A^^, and activation 
energies, are given in Table 2 (in column a) . These have been 

taken from Refs. 11 and 14 and are improved upon by using an 
optimal parameter estimation technique. 
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Table 1 


Kinetic Scheme in Terms of Functional Groups 


Polycondensation (Amidization) : 


A + B 


Imidization (Curing) ; 


Hydrolysis of Polyamic Acid (Degradation) : 


A' + B 


4) Hydrolysis of Anhydride Group : 


A + W 


5) Transformation 


Table 1 contd. 
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(Table 1 contd. . .b) 
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— OH 

— OH 


O 

11 

C OI I 

C — NH 

11 L 
o 



Table 2 


Kinetic Parameters 


Rate Constants 

* 

E. 

1 

( J/mol) 

* @ 

A. 

1 

(Initial Guess) 

(b) 

A. 

i 

(Best-fit 

Values) 

ki 

3 

(m /mol-s) 

2 . 93x10^^ 

2 @ 

7.42x10 

1.16x10® 

^1 

(s"^) 

1.03x10^ 

1 . 12x10® 

1 . 12x10® 

^2 

(s“^) 

1.09x10^ 

3.97x10® 

3 . 97x10® 


(m^/mol-s) 

1.50x10^ 

2 . 08x10"'^ 

2 . 08xl0“® 

k' 

(m^ /mol-s) 

1.03x10^ 

2.40x10^ 

2.40x10® 


(m /mol-s) 

5.27x10^ 

5.86x10"^ ® 

1.03x10® 

^4 

(s~^) 

7.62x10^ 

8.65x10® 

8 . 65x10® 

^5 

(s~^) 

9.63x10^ 

2.60x10"^ 

2 . 60x10”^ 

k' 

5 

(s'^) 

9.63x10^^ 

2.60x10"^ 

2 . 60x10"^ 


* Ref. 14 
@ Ref. 11 
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Before proceeding further, all the molecular species present 
in the reaction mass must be identified. The kinetic scheme of 
Table 1 suggests six sets of such species and these are defined in 
Table 3. In this table, -I- represents any internal group (-AB-, 
-A B- or -A B- , see Table 1). Each molecular species, e.g., 
^n AA' this table clubs together a host of molecules having the 
same number of repeat units, n, but having different chemical 
structures at the internal linkages. Thus, the concentration, 
[Pn aA^ ' represents a sum of the concentrations of all 
molecular species, P . having n AA units and n-1 BB units 
(with both the ends being an unreacted dianhydride) , but differing 
from one another in the sequential arrangements of internal units 
from among -AB-, -A B- and -A^B-. These molecules range from one 
having all the 2 (n-1) internal groups of the -AB- kind, to one 

ic 

having all 2 (n-l) of the -A B~ kind, etc. Such a clubbing is 
necessary to keep the model tractable. Thus, the concentration of 
P for example, is given by 


^^n,AA^ 


~ I '•^n,AA, i^ 

all species,! 


, etc , 


( 2 ) 


t b 

The k moments (k = 0,1,2,...) of the various molecular 

species are also defined in Table 3. Again, the summation over all 

species having the same length, n, b\it different internal 

* # . 

distributions of the -AB-, -A B- and -A B- groups, is accounted 
for by the use of the clubbed species (see Eq. 1) . The number and 
weight average molecular weights, and M^, for this complex 

system have been obtained starting from the fundamental 

8 9 

equations ' 
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Table 3 


Definition of Various Molecular Species, Moments 


and Average Molecular Weights 


* 

Molecular Species 


Symbol 

(a) 

Representation 

(schematic) 

Number of 

7\A 

BB 

AA' 

_I_ 

P 

n, AA 

AABBAA 

— ^AABBAA 

n 

n-l 

0 

2 (n-l) 

’^n,BB 

BBAABB 

— BBAABB 

n-l 

n 

0 

2 (n-l) 

P 

AABBAA 

BBAABB 

n 

n 

0 

2n-l 

n, AB 







p 

A'ABBAA 

^AABBAA 

n-l 

n-l 

1 

2 (n-l) 

n, A' A 







P 

A' ABBAA 

^AABBAA' 

n-2 

n-l 

2 

2 (n-l) 

n, A' A' 








A'ABBAA 

BBAABB 

n-l 

n 

1 

2n-l 

n. A' B 








* The total number of AA, AA' and BB (both at ends and in the 
inside) units in a molecule are listed above. The internal groups 
-I- (indicated as AB in column a above), could be either -AB-, 
-A*B-, or -a'^B-, as defined in Table 1. Thus, for example, 
clubs together several species having the same length and end 
groups, but with different internal groups. 

Table 3 contd. 
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(Table 3 contd. . .b) 


Moments of Various Molecular Species** 


k 


n=l 


n, AA 


BB 


AB 


. y n’" P 
L n. 


n=l 


. y n'^ P 

L n, 


BB 


AB 


n=l 


^A'A 


1 ^ ^ ,, 

L n, A' 

n=l 


A 


A'A' 


A'B 


I '■n.A'A' 


n=l 

00 


I ^n,A'B ; k - 0,1,2, 


n=l 


** P etc., indicate concentrations here. Again, the moments 

n f x\j\ 

club together molecules with different internal groups ( -AB- , 

"ic it 

-A B- and -A^B- ) . Note that k is a superscript and not an' 

exponent , 

Table 3 contd. 
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(Table 3 contd....c) 
Average Molecular Weights 


M 


n 


[( 


w 


^ux3 


AA ' "BB'^'^AA " '^BB ' “AB ' '^A'A ' '^A'A' 


I a: 


I A, 


•( 




^AA^^BB ^A'A ^ ^A'B ^ ^^A'A'^ 

^ ^A'A'^ + '^A'A^^A'A ^A'B ^ ^^A'A'^^ 


0 

''BB '■ '‘■AA 


w„^{A“, + A 


A'B’ 

0 

A' A 


r-)iO . -i ^ 'iO 

^T->n A^^ -h At 


BB 


AB 


+ A 


"a^r] 


r 2 r .T 2 ^2 

L ^AA^^AA ^BB 

0 


^^BB "■ ^BB^ '"bB^^AA ■ '^BB 


+ a: 


^ ^Ia} " 2 w^Wgg{A^ H- A^^ - A^ - A^g} 

+ 2 ^""aA ■‘' '^BB^ ^^A'A^^A'A ~ ^A'A^ 

^AA^ ^^A'A' " ^A'A'^ ■*■ ^^A'A ^BB^ ^^A'B ~ ^A'B^^ 


A' A' 


^^AA ■‘' '‘^BB^ ^^AB ^A'A ^A'A' ^A'B ^^A'A 

„.,1 -.0 ,.0 ,.0 1 2 ,0 
- 2 A^,^, - 2 A^,g + A^,^ + A_^,^, + A^,^} + 


A'BJ 


A' A A' A 


M 


w 


/o \ 2^0 / \ 2 ^ 0-1 

+ (2w-,- - w--) A-,,, + (w + w ) A-,J 


[('^AA ■*■ '^BB^ ^^AA ^BB ^AB ^A'A ^A'A' ■‘' ^A'B^ 

0 


'^AA^^BB * ^a'a ■*■ ^a'R ^^a'a'^ " '^rr^^ 


A' A 


A'B 


A'A' 


'BB^'^AA ^ ^A'A 


■vO 1 f ^0 A® . O-^^ \1 


A' A ^ A'A 


A'B 


A'A' 


PDI 


M 


W 


M.. 


n 


X 
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X 


M 


n 


= I 


(number fraction of species) 


all species 


(molecular weight of species) 


M 


w 


- ^ (mass fraction of species) x 


all species 


(molecular weight of i'' species) 


(a) 

(b) 
(3) 


The final expressions are given in Table 3. 

The balance equations for a well -mixed batch reactor for the 

various groups shown in Table 1 are now developed. These can 

easily be written and are given in Table 4. The concentration of 

the internal linkages, [-I-] , is the sum of the three 

concentrations, [-AB-] , [-A*B-] and [-A^B-] , and the balance 

equation for [-I-] can be obtained by summing up the corresponding 

three equations. The simplified equation for [-I-] is also given 

in Table 4. These equations can be integrated for isothermal 

operation (for a given set of rate constants and initial 

conditions) to obtain the concentrations of any group at any 

moment of time, t. The overall fractions of -AB- and -A B- in -I- 

are defined as f^ and respectively. These can be computed 

at any value of t, and the fraction of -I- which are of the kind, 
if 

-A"B-, can be obtained using 


fit = l-(f +f*) 
^A^B AB A B^ 


(4) 


The balance equations in Table 4 can give only partial 
information about the system. Mass balance equations for the 
individual (clubbed) species described in Table 3 need to be 
written before the equations for the moments of the chain length 
distribution can be obtained. These are summarized in Table 5. In 
these equations, great care has been taken to account for all the 
possibilities of formation and depletion of the individual 



Table 4 


Mass Balance Equations for the Groups in Table 1 


d [~A] 

= -k [-A] [-B] + k' [-AB-] - k. [-A] [W] + k;[-A'] 

dt -L 4 4 

d [-B] 

= -k [-A] [-B] + k' [-AB-] + k;rw] [-AB-] 

dt i 1 3 

d [-A'] 


= k' [W] [-AB-] + k. [W] [-A] - k;[-A'] 

dt ^ 4 4 

d [-AB-] ^ 

= k [-A] [-B] - k' [-AB-] - k„[-AB-] + k' [W] [-A B-] 

dt 12 2 

- k^[W][-AB-] - kg[-AB-] + k^[-A^B-] 

d [-A*B-] ^ 

= k2[-AB-] -k'[W][-AB-] 

dt 

d [-A^B-] „ 

= k^I-AB-] -k'L-A’^B-] 

dt ^ ^ 

d [W] 

= k2[-AB-] -k^[W][-AB-] -k^[W][-AB-] - k^ [W] [-A] 

dt 

+ k^[-A'] 

Some Other Quantities 


d [-I-] 
dt 


{ [-AB-] + [-A*B-] + [-A^B-]} 

dt 

k^[-A][-B] - k^[-AB-] - k^[W][-AB-] 

Table 4 contd. 


23 



(Table 4 contd 


b) 


'AB 


f ^ 

A B 


^A^B 


[-AB-] 


[-I-] 

[-A*B-] 

[-I-] 

[-A^B-] 

= 1 - f - f * 

r T 1 AB ^A B 
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Table 5 


Mass Balance Equations for the Molecular Species 

in a Batch Reactor 


dP 


n-1 


n, AA 


dt 


^^1^1 I ^r,AA^n-r,AB 


p OD 00 

^^l^n,AA[ Z ^m,AB ^ Z 


AB ' L “m,A'B 


r=l 


m=l 


m=l 


+ 2 


00 ^ 
Z ^it\,BB J 


^l^AB 


p 00 00 

_ Z ^m,A'A ^ ^m,AA 


H 

If 

e 

m=n+l 

■' f (n-1) P ' - 

■1 AB^ ^ n,AA 

2k' f W 
3 AB 


+ 2> P ,, + 

m=n+l tn=n 


CO 1 

Z ^m,ABj 


n,AA ^^4’™'^n,AA 


^4^n,A'A 


dP 


n-1 


n,BB 


dt 


^^1^1 Z ^r,BB^n-r, 


AB 


r=l 


00 -1 r “ . 

y p + k'f,„ 2) p„ 

L m,A'Aj 1 AB|_ L m, 


BB 


m=l 


m=n+l 


00 


00 


^m, AA 

^ z p 

. 

m= 

= 1 

m=l 

00 


00 

V P 

_ _ 4 * 

y p 

L m, AB 

L m, 

m=n 


m=n 


p 00 00 

2k'f.„(n-l)P„ 2^ P^^gg + X Pm,AB 


"I AB 


n,BB 3 AB 


m=n+l 


m=n 


I ^m,A'B] 


2k'£^W(n-l)P^^BB 


m=n 


Table 5 contd. 
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(Table 5 contd. . . .b) 


n, AB 


= 4k. y P .^P 
lA, m,AA n- 


m+l,BB ^1^1 I ^m,AB^n-m,AB 


]r p 

1 n,AB 


[ ra 00 

^ ^m,AA ^ 

m=l m=l 


m,BB ^m. 


L m,AB L m,A'A 
m=l m=l 


UJ 

X ^m,A'. 


00 

I ^tn,A'Bj 


H^AB I H- 2P^_gg + 2P^^^ ^ 


1 AB Z. ’• m,AA m,BB 

m=n+l 


m=n+l 


’*' Zi'aJ ^ k-ofaTa^^^n 1) P 


k,WP 


m,AB ' ^m,A'A^ -3"AB”^"“ "^"n,AB ^4'’"n,AB 


■^^4^n,A'B 


n,A' A 




CXJ 00 

,A'A_ ^m,BB ^ I ^m,AB 


La m, BB La m. / AB Lu 

m=l m=l m=l 


CO -. 

I ^m,A'Bj 


+ 2k£A^^ ^m,A'B^n-m,AA ^l^ll ^m, A' A^n-m, AB 


^^i^AB ^n,A'A ^1 ^Ab[ I ^m,A'A ^m,A'A' 

m=n+l ni=n+l 


I ^m.A'B] - 2ki*AB<''-l>" 

m=n 

00 00 

Z ^m,AB ^ Z ^m,A'AJ ~ ^ 


P 00 

^^3^AB ^”'"^^^^n,A'A ^3^AB^L ^Z ^m,AA 


m=n+l 


^WP + 2k. WP 

4 n,A'A 4 n,AA 


m=n m=n+l 

^4[ ^^n,A'A'" ^n,A'A] 


Table 5 contd. 
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(Table 5 contd. . . .c) 


n, A'A' 


k_A. y P 

1 1 A m,A' 


A^n-Tn,A'B ^^^n,A'A' 


jr 00 

^n,A'A' ^3^AB^[ Z ^m,A' 


I ^m.A'A Pm,A'A'] +'"4"^n,A'A - 2k4Pn,A,A, 


m=n+l 


m=n+l 


n,A'B 


[ 00 00 . 00 -1 
^m,AA I \,AB Z ^m,A'Aj 


^l^lZ ^m,A'B^n-m,AB ^m,A'A^n- 


m+1 , BB 


HfAB<2"-l>A,A'B * 


P 00 00 

^1^Ab[ I ^m,A'A I ^m,A': 


m=n+l 


m=n+l 


^ ^m,A'A'] 




m=n+l 




[ 00 00 00 

Z ^m,AB ^Z ^m,BB Z ^m,A'A 
m=n+l m=n+l m=n+l 


'I ^m,A'B " 2 


00 -I 

Z ^m,A'A'_ 


+ k^WP - k;P 

4 n,AB 4 n,A'B 


m=n+l 


m=n+l 
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(Table 5 

dW 

dt 


contd . . . . d) 


p 00 

- ’^3«Ab"[I " ^b.A'b) 


4 - P 4 - P 4 - P \ 

AA n,BB ^n,A'A n,A'A''^^ 


n=l 


k^w 


p 00 

[ I 


AA 


n=l 


+ P + p 

n,AB n 


,A'A^] 


+ k; 


I ^’^n,A'A ' "n,A'B 


+ P 


n=l 


+ 2P 


n 


,A'A' 


+ (k f - k'f * W) 
^ 2^AB 2^A B ^ 


[ I (2n-2){P„_ 


AA 


n=l 


n=l 


00 

^n,BB ■*■ ^n,A'A ^n,A'A' ^ Z ^^n,AB 

n, A'B^] 


+ P 


a. = 0 if n = 1 
= 1 if n > 1 

b. f^g and fj^*g to be obtained from integration of (all) the 
equations in Table 4. 

c. notation [,] for concentrations are not being used in this 
table for the sake of brevity. 
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molecular species which the kinetic scheme (in terms of groups) in 

Table 1 represents. The terms in the first equation in Table 5 

(for are explained to illustrate the several 

assumptions and concepts involved. The first term, 

n-1 

^r,AA^n-r AB' ^^P^esents the formation of t)y the 

reaction between A and B end groups of the shorter P and 

^n-r,AB i^olecules, the factor of 2 accounting for the two 
possibilities (2 A groups of ^ reacting with the single B 

group of • The next three terms involve the depletion of 

^n AA reaction of two of its end A groups with the end B 

groups of all P^ P^ and P^ gg molecules. The reaction 

between P with P _ involves four possibilities, which is 

li / jc\jc\ in f JlSJis 

reflected by the multiplication factor. The next three terms 

involving k^f^ account for the formation of P^ ^ by the scission 

(due to the reverse Rxn. 1) at an appropriately placed internal 

-AB- link of a higher length P 7 v, 7 v, P ,, or P molecule (m > 

n or n+1) . While only a single molecule of P is generated when 

P and P undergo this reaction (the other molecule being 

m f jc\, jc\. m f 

P and P respectively), there are two Ideations of 

m-n,A'B m-n,AB' j < 

-AB- in P which could lead to the formation of P (the 

m, AA n, AA 

other molecule produced being - The factor, f^g, used in 

these terms is the probability that the appropriate internal 
location at which this chain scission is to take place, is indeed 
an -AB- linkage rather than an -A B- or -A^B- . It is being assumed 
that the overall fraction of -I- groups being -AB- is the same as 


the local probability of finding an -AB- group at the right 

location in any of the larger molecules a^A^ ^m AA AB^ ' 

Such an approximation that the overall probability is the same as 


W' 



the local probability, independent of the length, m, is quite 
commonly made in polymer reaction engineering, ^ particularly for 
the polymerization of phenol and melamine formaldehyde.^^ The next 
term, -2k£f^g (n-1) ^ , represents the depletion of by the 

reverse step of Rxn. 1 (Table 1) at any of its 2 (n-1) internal 
locations which happen to be of the -AB- kind. Again, the 
probability that any of these internal locations is -AB-, as 

■*. J1 

opposed to -A B- or -A^B-, is assumed to be f^, the overall 

fraction of internal locations which are of the -AB- type. The 

term, -2k2f^gW (n-1) P^ can be similarly explained for Rxn. 3. 

The last two terms in this equation correspond to the depletion 

(two possibilities at the two -A end groups) and formation (only a 

single possibility) of P ,, and P by the forward and reverse 

steps of Rxn. 4 in Table 1. All other equations in Table 5 can be 

similarly explained, using the above discussion. The equations in 

Table 5 were checked out by two independent derivations. In 

addition, using k. = k^ =0 (i^l) , it was found that these 

11 

equations reduce to those for simple AA + BB reversible 
polymerization (without formation of the condensation byproduct) . 
These give us confidence on the correctness of these relatively 
complex equations. 

The equations in Table 5 can be appropriately summed up 
(using the definitions in Table 3) to give ordinary differential 
equations (ODEs) for the several moments (k=0, 1 and 2) . Several 
summation identities given in Chap. 1 in Ref. 9 (besides a few 
others which had to be developed) proved helpful. The final 
equations for the moments are given in Table 6. It is observed 
that the equations for several of the second moments require third 
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Table 6 


Moment Equations in a Batch Reactor 


dA. 


AA 


dt 


+ X^,^- X°,^) 
~ 2)c^WX^^ + k^X^, 


4 AA 
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,0.0 > -.1 .1 .^1 ^0 

^AB^A'A^ ^l^AB^^^AA ^^BB ^A'A ^A'B ” ^AB 


- 2A^ - 2A^ 

AA BB 


4'A- 4’B> " >'3*AB”<2^iA + ^I'A 


^ 1 0 .0 > 1 TT^ 0 1 , , 0 

- ^AB - - VJ'ab + ''4^A'B 


dX 


A' A 


dt 


^1 ^^^AA^A'B 


2A° A° 
^^BB^A' A 




-1 .1 .0 

^A'B ' ^A'A ^A'A 




^A'A ^A'A ^^AA^ ^4^^^^AA ^A'A^ ^4^^^A'A' 

- ^A'A^ 


dX 


A'A' 


dt 


^l^A'A^A'B " ^^i^AB^^A'A' ' ^A'A'^ ^3^AB^^^A'A 

^A'B “ ^A'A^ ^4^^A'A " ^^4^A'A' 
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(Table 6 contd. . .b) 

0 ■ 


dX 


A'B 


dt 


0 .0 


“ ^^L.^A'B " ^A'A^A'B^ M^AB^^^A'A' 


0 .0 


+ X 


A' A 


A'B 

1 




^^A'A'^ ^B^AB^^^^BB ^AB 


+ + 2X^,^, - 2X°g - X^ - X^,^ - 2X°,^, - X°,g) 


4 A'B 


dX 


AA 


dt 


^1 ^^^AA^AB 


4xLx° 


AA'^BB ^^AA^A'B^ ^1^AB^°'^^AB 


+ 0.5X 


A'A 


°-5^iB - °-^4 'a) 


^^B^^AB ^^AA 


^AA^ ~ ^^ 4 '^^AA ^4^A'A 


dX 


BB 


dt 


“^^AA^BB " ^^A'A^Bb' ■*■ ^AB^ ^ ^ ^AB 


* °'^^A'B ^BB * ^BB * * °'^^A'B* 


* k'f^W(0.5xfg . 0 . 5 X ^,3 


4b * 


X 0.5X-,g) 


dX 


AB 


dt 


- k f4X° X^ + 4X° 4A° 2>° >1 

^l^^^AA^BB ^^BB^AA ^^AA^BB ^A'A^AB ^^AA^AB 


^A'B^AB 

^ 0-5 4'A 


0 1 > 
'■BB AB^ 


k' f 

(X^ 

1 AB 

^ AA 

.1 

.1 

^AA ~ 

^BB 


2 

'■BB 


AB 


- O.SX^,^) -k^WX^ ■’■ *^4^A'B 


0.5X1^ -0.5X1^) 
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(Table 6 contd. . . .c) 


1 AA A'B ^A'B^AA ~ ^^BB^A'A ^A'B^A'A 

- ^i-A- * O-S^A'B ♦ 1-5 ^A'a' + '^3^Ab”'^L 

+ 0.5X^ - 1.5X^,^ - + 0.5X^ X I-S^'a' 

+ ^4^^^^AA ' ^A'a' '*' ^4^^^A'A' "-^A'A^ 


0 . 1 


0 ,1 


^^A' A' 


‘^^A' B 


’^l^^A'A^A'B ■*■ ^A'B^A'A^ ” ^’^l^AB ^ ^A' A' " ^A'A' 

. k'f^^W(0.5A2,^ , O.sxl^ - Xl:^, . Xl^, 

- O.sxl^ . ^ V4'A - 2J-4^i'A' 


k (2X^ + 2X^ X^ - X^ 

^l^^^A'A'^BB "^^BB^A'A ^^BB^A'A ^^AA^A'B 

\0 \1 ^0 \^\ lr/-P f n TCTk^ 

^A'A^A'B ^A'B^AB' ^1^AB^°'^^A'A ^A'A' ^'^^A'B 

X 0. 5X^,3 - O.SX^,^ - xj,^,) X k^f^W(X^3 X 0 . SX^ 

X O.SX",^ - xj,3 X X^,^, - X^3 - 0.5X^,^ - 0.5X^ 

V'A'* * ’'4^'aB " ’^4V'B 


, ,-^0-2 x\l-,l A-^0 .2 

^l^^^AA^AB '^^AA^AB ‘^^BB^AA 


2X° ) 

^^A'B^AA^ 


, , . ,1-3 ^1.3 _ 4 - 3 ^ - 2 2 

’^l^AB^ ^AB 3 ^A'A ■ 3 ^AA ^AA °-^^AB 


,4 r-.2 1 ^1 1-11,1^ 
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2k^WX^ X k4x2,^ 
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(Tabie G coiiLd. . . .d) 


^1^AB^3 ^AB "'' 3 ^A'B " 3 ^BB "'' ^BB ''' °*^^AB 

+ 0.5X^,g + I ^BB + ? ^AB I ^A'B^ ■" ^S^AB^^I ^AB 

'*' 3 ^A'B 3 ^BB ^BB ■*' °'^^AB ^’^^A'B ■’' 3 ^BB 

^ I ^ab "■ I 4'B^ 


^l^^^AA^BB ^ ^^BB^AA ®^AA^BB ^ ^^AA^BB 




^^Ia4b 


2X° X^ - 2X° 0 2 
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A° >2 „ 1 1 , 2,3 ^1 ,3 
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1 .3 


+ f K'n - - °-^K.u * k ^In *i 


3 A'B '"BB 


A'B 3 BB 6 A' A 


1 -1 ^ 2 ,3 4-3 -2 1-1 1 a ^ 

6 A'B 3 ^AA ■ 3 ^AB " ^AA 3 ^AB ^ 3 ^AA‘ 

,,,. ut,2 .3 5.3 1-3 a nc^2 

^3^AB^^3 ^AA 3 ^AB 3 ^A'A ^AA ° ’ ^^AB 

0.5Xa,a + 2 ^aA 6 ^AB ^ 6 ^A'A^ ~ ^4^^AB 

+ ^4^A'B 


<'A 


1 _a 12 ^.^1^1 ^-0 2 
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, -j > ^ ^ 9 A ^ ) + k ' f (— X^ 
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- X^ + 0.5X^ + 1.5X^,^ + 3 X^ + -g- X^ + -g- 

+ k4W(2X^ - + ^b^^A'A' ” ^A'A^ 
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(Table 6 t;ont:d. . . .<') 


A' A' 


0 , 2 
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1 . 1 
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moments. Thus, we have a hierarchy of infinite equations. A 

standard technique in polymer reaction engineering^'^ is to break 
this hieiaichy by using closure conditions. The closure equation 
which has been found to be quite successful for nylon 6, 

polyester, polyethylene and polyphenylene oxide^"^ reaction 
systems, is tried here. This is also included in Table 6. 

The following additional checks were made on the equations 
dsrived, to ensure that these were free of errors: 

1) In Table 7, several groups (end as well as internal) have 

been written in terms of the various moments of the molecular 
species. Using these equalities with the equations of Table 
6 , we can obtain expressions for the rate equations 

analytically . We obtained the same equations as given in 
Table 4 (which were obtained using a functional group 
description, rather than a molecular species approach) . 

2) In a batch reactor, the total number of repeat units, AA, BB 
or AA' , should be the same at any t as at t=0. Using Table 3, 
it is possible to obtain an equation for the total number of 
each of these repeat units in terms of the various moments. 
It is found that the rate of formation of each of these 
repeat unit is zero, as expected. 

These two checks were sufficient to confirm the correctness of the 

mass balance and moment equations in Tables 5 and 6. 

The equations in Tables 4 and 6 form a complete set of ODEs 

(initial value problems -- IVPs) , which can be integrated for any 

given set of initial conditions (feed values) . The numerical 

2 8 

technique used to integrate these equations is Gear s algorithm. 
The NAG library routine, D02EBF, which has a built-in step size 
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control algorithm and is particularly useful for stiff systems, is 
used (with a tolerance, TOL, of 10 for this purpose on an HP 
O000/8t0 .q computer system. Changing the value of the parameteV, 
TOL, led to almost identical results. The two checks described 
above were also made with the computer code to ensure the 
correctness of the program (referred to as the simulation 
program) . 

Our set of modeling equations and our approach differ from 

the two models currently being used. In one of these, the 

stochastic theory of chain branching processes of Gordon et 
2 9 3 0 

al. ■ ' has been used. It is difficult to accomodate mass and 

heat transfer aspects in such theories. The other approach uses 

23 24 

the results of Kilkson, ' who develop expressions for only the 
final values of the average molecular weights for irreversible AA 
+ BB condensation polymerization, in which all of the BB monomer 
is charged into a reactor, and AA is then added at some rate till 
the desired overall stoichiometric ratio is attained. It is 
inappropriate to use these results for the more general, 
reversible kinetic scheme of Table 1, particularly to stage 2 
polymerization where reversibility effects and side reactions are 
important. 

Sensitivity studies were made using the computer program for 
the polymerization conditions: 

[-A] Q = [-B] Q = 400 mol m'^ 

I-A'Iq = [-AB-]q = [-A*B-]q = [-A^B-]q = 0 mol m'^ 


[Wlo 

= 

55 mol m 



[ = 200 mol m’^ 

0 " 

^ ^AB^ 

0 

^ ^A'aL ~ *• ^A'A'L 


; k =0,1,2 
^ ^A'B^ 0 "" 


0 mo l m 


-3 


; k = 0,1, 
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The several rate constants (not the Arrhenius parameters) were 
varied by ±3 0 %, one at a time to study their effects. It was 

found that the results are sensitive to only two parameters, 
and k^ . A two-parameter curve-fit was then performed [keeping all 
other rate constants at their guess values (Table 2)], to obtain 
best -fit values of k^ and k^ at 27 *^C. In order to carry out this 
'tuning' study, the computer code for simulation was combined with 
a nonlinear parameter estimation program. The Box^^ complex method 
was adopted for this purpose. 

The Box complex procedure obtains values of the parameters, p 

( = [p^ / P 2 ' • • - ' Pq] ^ ' which minimize some objective function, E(p) , 

using a patterned-search technique. The objective function is 

chosen as a weighted sum-of -square errors between values predicted 

by the model and some (scarce) experimental data, reported in the 

11 

open literature. It is given by 


E (p) = E (k^, k^) = ^ 

i = l 


N , M . - M . ^2 

n,l,expt n,l, theor 


M 


n , 1 , exp t 


N r M . 

^ w . w , 1 , c X p t 


- M . .2 

W , X , tlieor- I ( 6 ) 


1=1 


W , 1 , ex p t 


In the above equation, the values of and are those at t=t^ 
(s36000 s) , the final reaction time (after which almost no further 
reaction takes place) , N is the number of experimental data points 
available for and and w^ is the weightage factor assigned to 
the i*^^^ data point. Superscripts, expt and theor, represent 
experimental and the corresponding theoretical values. The tuned 
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rate constants are then used to generate frequency factors for 

and , assuming that the activation en^^rgies are the same as the 

guess values in Table 2 . These values are given in column b in 

Table 2. Using these parameters, results other than for M and M 

n w 

(required for parameter estimation, for which experimental data is 
available) at 27 °C are generated. The CPU time taken on a 
mainframe HP9000/850S supermini computer for the parameter 
estimation (18 iteration) was 90 s. 
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CHAPTER 3 


RESULTS AND DISCUSSION 

The Box complex program is used as described above to obtain 
the best-fit values of and at 27°C. The computational 

parameters used in this program are given in Table 8. Figures 1-4 
show some of the results generated using the Box complex code for 
parameter estimation {for different values of [P^ ^ + [P^ ^ 

(equimolar amounts) but with other conditions the same as in Eq. 
5}. Figure 1 shows the final values (at t^ = 36000 s) of the 
molecular weights, M and M , predicted by the model, using the 
tuned parameters (Table 2, column b) . These are found to compare 
well with the experimental data points of Kolegov et al. for 
different values of the (sum of) initial monomer concentrations, 
[P^ AA^ 0 '■^1 BB^ 0 ■ 2 shows how the error, E (see Eq.6), 

decreases with iteration number, while Figures 3 and 4 show how 
the two parameters, k^ and k^ , evolve with iteration number to 
their (near) optimal values. Figure 1 also shows that the optimal 
parameters fit the data better than their initial guesses, which 
were fairly good. It would have been better if actual experimental 
points were available as a function of time for each of the values 

^^1 AA^ o’*" ^^1 BB^ 0 ■ Unfortunately, only model predictions of 

11 

these properties have been presented by Kolegov et al. 

After obtaining the best-fit values of the parameters, p [ s 
k , k ] , these were incorporated in the simulation program and 
additional results were generated for the feed conditions of Eq. 
5, i.e., for [Pi^aa^o = ^^1,BB] " 200 mol m“^ . Figure 5 shows 

the build-up of the number and weight average molecular weights 
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Table 8 


Data Used in Optimal Parameter Estimation 


Weightage Factors in Eg. 6 



=16 

D 

= 32 

D 

10 

Values of Parameters of Box Complex Program 

N =2 
a =1.3 
(3 = 10 "^ 

r =5 
5 = lO"^ 

Initial Value of Parameter, p 

p H [k^, k^] = [6.0, 0.04] 

Bounds on p 

0.0 :£ k ^ ^ 10.0 
0 . 0 k^ :£ 1 . 0 
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their initial values {Table 2 column a) . Experimental data 
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Fig. 2 Variation of the error vs. iteration number. 
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Fig. 3 Variation of k- with iteration number. 
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Fig. 4 Variation of k. with iteration number. 
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with time. These curves show the typical characteristics of a 
reversible polycondensation reaction in which the molecular 
weights reach an equilibrium value after some time. This kind of 
asymptotic behavior was also observed for other values of 

studied in Fig. 1. Similar qualitative 
behaviors have been predicted by Kolegov et al.^^ under similar 
conditions. 

The effect of the two rate constants, k, and k„ , on M and M 
are shown in Figures 6 and 7. These figures show clearly that 
and are quite sensitive to these two parameters, but that the 
polydispersity index, PDI, is very close to 2 . 0 , as shown for the 
reference conditions (optimal values of and and feed 

conditions of Eq. 5) in Fig. 8 (the effect of k^ and k^ on PDI is 
minimal) . Other parameters have been found to have no significant 
effect on and M^. Increasing k^ while keeping k£ unchanged 

speeds up the first reaction as expected, and also leads to higher 
molecular weight product at equilibrium due to increase in the 
value of its equilibrium constant. A more interesting observation 
is that the molecular weights decrease as k^ is increased. This is 
because the increased hydrolysis of the anhydride groups at the 
ends of the chains to form -A' associated with higher k^, leads to 
a depletion of -A groups and this suppresses the chain lengthening 
process associated with the first reaction (-A' does not 
participate in any chain lengthening reaction since the third 
reaction is irreversible, as shown in Table 1) . The lack of 
sensitivity of results to the rate constants associated with the 
other reactions is because they are relatively unimportant at 
27°C, the usual temperatures of polymerization in the first stage 
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Pig. 6 Effect of varying k on M (dotted) and M (solid) . 
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reactor being studied. It is expected that at higher temperatures 
(around 130 C typically encountered in stage 2 reactors) , the 
other reactions would play a more important role, and the results 
may be sensitive to variations in the associated rate constants. 
It is clear that these rate constants should be 'tuned' against 
experimental data taken on second stage reactors whenever they 
become available in the open literature. Till that time the 
general mathematical framework developed here remains un-tuned for 
the second stage reactors, and can be used only to 
provide qualitative trends (as done later in this paper) . 

Figure 9 shows the conversion of the monomer, P , with 

JL f XTuC^ 

time under reference conditions. It is found that the conversion 
of the other monomer, is almost indistinguishable from that 

of P^ indicating that the fourth reaction in Table 1 (which 

depletes -A groups but not -B groups) does not contribute 
significantly to monomer consumption in the first 4 seconds in 
which most of the monomers are consumed. The rapid depletion of 
monomer is typical of step growth polymerizations. It may be 
emphasized that the fourth reaction does, indeed, affect the 
molecular weight at later times (this explains the decrease in 

and M as k increases, in Figure 7) , which is associated with an 

W Tt 

increase in [-A'] with t (see Figure 10). The concentration of -A 
functional groups decreases at a slightly slower rate than 

'■^1,AA^ ■ 

The average molecular weights of PAA can be controlled by 
varying the stoichiometric ratio, R = ^ ^ small 

vax'iation in R away from unity results in a significant drop in 
the molecular weight, as shown in Figure 11. This behavior is 
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similar to that exhibited by (irreversible) AA + BB step growth 
polymerizations, S which is not surprising since Rxn. 1 in Table 1 
is Che mam reaction under reference conditions. Again, the PDI 
attains a value of 2.0 quite early. The reaction (in stage 1) is 
also quite sensitive to the initial water concentration.^^ Figure 
11 shows how the average molecular weight changes when [W] is 
increased from the reference value of 55 to 66.0 mol (1.2 
times the reference value) . The importance of proper metering of 
the monomers and control of their purity is brought out by this 
diagram. Figure 12 shows how the fraction of imide links and of 
the unfavoured conformer, f^*^ and f^#^ respectively, vary with 
time. Both these fractions are negligible at reference (first 
stage) conditions, and f^ rises very rapidly to a value near 
unity. This confirms that Rxns. 2 anc| 5 are relatively unimportant 
at low temperatures. The variation of the acid group 
concentration, [-COOH] , with time is shown in Figure 13. These can 
be experimentally measured by titration. ' The increase in [-COOH1 
with time compares qualitatively with the experimental results of 
Johnson on a different system (his equivalent weights are 
evaluated such that they show the opposite trends than shown by 
C-COOH] ) . 

The moment closure approximations (Table 6) were changed by 
incorporating a multiplying factor on the right hand side of the 
equations. No significant differences were found in all the 
results at 27 °C when a factor of 1.2 (or larger) was used. This 
insensitivity of results justifies the use of these closure 
conditions and is probably because under these conditions, only 
the first reaction in Table 1 contributes significantly (and the 
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(dotted), with time for the reference conditions. 




Fig. 13 Variation of [-COOH] with time for the reference 




the reference conditions (but with T different than in 










closure conditions have been well tested for this reaction^) 

Figures 14-16 show how the values of [“A*B-] , [-AB-] , [-COOH] 

and at large values of t (= 3 600 0 s) , vary with temperature. 

The initial conditions are the same as given in Eq. 5. Again, 

7 . 

Johnson has provided some experimental results similar to these 
(for PAA prepared from a different diamine) . Figure 14 shows that 
much higher concentrations of imide links are attained at higher 
temperatures, at the cost of the amide links {-AB-) . Such a 
phenomenon has been predicted qualitatively by Johnson , ^ and 
reflects the importance of Rxn. 2 in Table 1 at high temperatures. 
Figure 15 shows that the concentration of acid groups decreases at 
high temperatures as they get . consumed by the imidization 
reaction. This is exactly opposite to the trend shown by the 

. 7 

equivalent weight vs. T plot of Johnson. 

Figure 16 shows how at high values of t (= 36000 s) varies 
with temperature. This diagram also shows how an increase (or 
decrease) in is associated with a simultaneous decrease (or 

increase) of the concentration of -A' . Figure 17 shows (at t = 
36000 s) vs. T for several combinations of the reactions of Table 

1. Curve a shows the results when only Rxn. 1 (Table 1) is 
considered. Higher temperatures lead to higher values of 
(relative to k^) , and lead to lower chain lengths for the 
exothermic Rxn. 1. Just by adding on Rxn. 4 to this reaction, we 
find (curve b) that the value of is reduced by almost three 
orders of magnitude. This is because Rxn. 4 depletes the -A end 
groups (converting them to -A') and thus reduces the chain growth 
process of Rxn. 1. It is obvious (see curves b and c) that Rxns. 

2, 3 and 5 do not contribute significantly till about 80 °C. The 
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decrease of followed by an increase and then, again, a 

decrease, as T increases from about 25°C to 90°C is, clearly, due 

to the effect of Rxn. 4 on Rxn. 1 (see curve b) . The reverse step 

of Rxn. 4 dominates in 60°C ^ 90*^C while the reverse step of 

Rxn. 1 dominates for T > 90 °C. In fact, the local minima/maxima 

is not observed if only k^, and k^ are nonzero (see curves a, b 

and d) . At T > 90 °C, Rxns . 2, 3 and 5 also affect M (see curves 

n 

b, c and e) . W produced by Rxn. 2 leads to chain scission through 
Rxn. 3, thus giving lower molecular weight product. It may be 

, 7 

mentioned that Johnson has found that the viscosity of polymer 
formed at large values of t decreases as the temperature 

increases. However, he does not study very high temperatures and 

so it is not too clear if the behavior shown in Figure 17 is 

indeed observed experimentally. If detailed experiments become 
available, one can 'tune' the rate constants associated with Rxns. 
2, 3 and 5 and get better predictions. 

Figures 18-20 show some results for simulating the second 
stage (curing) operation (earlier results were for the first stage 
at different T) . The initial conditions for generating these 

results are the same as at the end of stage l (t = 36000 s in 
stage 1, at 27°C) ; 

Stage 2 
at t = 0 : 


[-Alo 

[-A'] 

[-BIq 


0.47 

X 

lO"'^ 

mol 

m 
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[-AB-]q = 

0.40 X lO"^^ 

mol 

-3 
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0.21 
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10"^^ 

mol 
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other values (moments) can be supplied on request. The temperature 
is selected as 130 C (isothermal) . It is assumed that the solvent 
is not removed at the end of stage 1 (i.e., a film is not cast and 
then cured) . Since little experimental data is available for the 
curing of PAA. films against which our results can be tested (or 
tuned) , and since the objective of generating these results is 
solely to test whether our model can, indeed, predict some 
important qualitative features expected in stage 2, removal of the 
solvent before curing was not effected. Indeed, several additional 
physico-chemical phenomena (e.g., removal of water through the 
solid film, dependence of rate constants on the 'viscosity' of the 
film, heat transfer across the solid film, etc.) need to be 
accounted for simultaneously in the present kinetic model before 
curing of PAA films at high temperatures can really be modeled. 

Figure 18 shows how and decrease with time (due to the 
chain scission reaction) .. However, these increase slightly after 
about t = 20000 s. This is associated with a simultaneous increase 
in -A' concentration, as shown in Figure 19. The continuous 
increase of f,*^ arid f,#^ with time for the same situation is 
shown in Figure 20. It is clear from Table 1 that the reverse step 
of Rxn. 4 starts playing an important role at high values of t, 
and converts -A' to -A (relatively speaking) , thus leading to 
higher and (through Rxn. 1) as well as lower [-A'] at large 
times. The importance of the reverse step of Rxn. 4 in Table 1, is 
to be noted. It has also been found that the molecular weights are 
relatively insensitive to variations of and [W] q in stage 
2 . 

The effect of introducing a multiplying factor of 1.2 on the 
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right hand sides of all the moment closure equations in Table 6, 
was also studied. The value of (at t = 3 6000 s) decreased by 
less than a few percent by using this factor in all the closure 
equations simultaneously. The difference was quite small if only 
one closure equation was modified at a time. This suggests that 
these closure equations can possibly be used for the second stage 
reactor. However, the sensitivity of results to the moment 
closure equations needs to be studied again after tuning the rate 
constants for Rxns .2,3 and 4 on experimental data for the curing 
process. 
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CHAPTER 4 


CONCLUSIONS 

A comprehensive kinetic model has been presented in this 
study. This model can easily be adapted to account for limitations 
of heat and mass transfer which are normally encountered in large 
scale reactors. A few of the rate constants have been curve-fitted 
against a limited amount of experimental data available at low 
temperatures, from the open literature. However, detailed 
experimental results, particularly at high temperatures, are 
required so as to get 'tuned' values for the other rate constants, 
in particular, Rxns. 2, 3 and 5 of Table 1. 
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